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The most simple and seemingly straightforward application of the photon blockade effect, in which 
the transport of one photon prevents the transport of others, would be to separate two incoming 
indistinguishable photons to different output ports. We show that time-energy uncertainty relations 
inherently prevent this ideal situation when the blockade is implemented by a two-level system. 
The fundamental nature of this limit is revealed in the fact that photon blockade in the strong 
coupling regime of cavity QED, resulting from the nonlinearity of the Jaynes-Cummings energy 
level structure, exhibits efficiency and temporal behavior identical to those of photon blockade in 
the bad cavity regime, where the underlying nonlinearity is that of the atom itself. We demonstrate 
that this limit can be exceeded, yet not avoided, by exploiting time-energy entanglement between 
the incident photons. Finally, we show how this limit can be circumvented completely by using 
a three-level atom coupled to a single-sided cavity, enabling an ideal and robust photon routing 
mechanism. 

PACS numbers: 42.50.Ct, 42.50.Pq, 42.50.Ex, 42.50. Dv 



I. INTRODUCTION 



Often described as the optical equivalent of Coulomb 
blockade JL> »2, , photon blockade [3j , in which the trans- 
port of only one photon through a nonlinear system is 
possible whereas excess photons are rejected, is one of the 
few photon-photon interactions that have been demon- 
strated experimentally to date [!4-*B1. Such nonlinear in- 
teractions at the single-photon level are inherently non- 
classical, and form the basis for quantum information 
processing with photonic qubits [7l410j. In particular, 
the demonstration of the photon blockade is considered 
a significant step towards the goal of single photon rout- 
ing '11' . Two mechanisms of photon blockade have been 
demonstrated in recent years. Both mechanisms rely on 
the coupling of an atom (or another two-level system such 
as a quantum dot [H]) to cavity-enhanced electromag- 
netic modes. In this work we establish that such an ap- 
paratus is inherently insufficient for the task of photon 
routing, and that a three-level system at least is needed to 
ensure a deterministic and efficient routing process. We 
demonstrate this by analyzing the simplest case of an 
input pulse containing exactly two photons, and deriv- 
ing the probability for a successful routing event, namely 
that the system will respond to one photon differently 
than to the other. 

The underlying conflict that limits the interaction of a 
two-level system with such a pulse, is the fact that the 
coupling to the electromagnetic mode dictates both the 
interaction bandwidth, and the memory time of the sys- 
tem. Thus, a pulse that is short enough to guarantee that 
the two photons arrive within the memory time of the 
system will have a bandwidth that exceeds the interac- 
tion bandwidth. Conversely, a pulse that is narrow-band 
enough to be included completely within the interaction 
bandwidth will be long enough to allow the system to 
respond to each photon as if it were the only one. 



The outline of the paper is as follows. In Sec. |TT] we 
present both analytical derivations and numerical cal- 
culations of the efficiency of two-photon routing based 
on photon blockade in cavity quantum electrodynamics 
(QED) [H H- The temporal behavior of the block- 
ade mechanism is essentially the same in the bad cav- 
ity regime as in the strong coupling regime, despite the 
fact that the underlying nonlinearity is different in the 
two regimes. We show that the routing efficiency is in- 
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herently limited to ~ 64% in both cases. In Sec. 
we study the possibility of exceeding this limit by using 
a time-energy entangled photon pair as the input pulse. 
This case seems especially relevant in light of the fact that 
such entanglement can be created by a two-level system 
coupled to a single electromagnetic mode (THl US], and 
so one could think of a two-stage process in which the 
first interaction creates the entanglement between the 
photons and the second interaction performs the rout- 
ing mechanism. However, our analysis shows that even 
ideal time-energy entanglement can increase the routing 
efficiency only up to ~ 77%. Time-energy entanglement 
generated by a two-stage interaction with a two-level sys- 
tem is even more limited and increases the efficiency only 
up to ~ 68%. Finally, in Sec. IV we study the config- 
uration of a three-level atom coupled to a single-sided 
cavity. Recent studies by Koshino et al. [17] and Gea- 
Banacloche et al. [TH] have shown that this configuration 
enables deterministic mapping of a photonic state to the 
atom. We utilize this scheme to construct an ideal pho- 
ton router that does not suffer from the inherent limits 
of photon blockade with a two-level system, and show 
that its efficiency can approach ^ 100% with realistic, 
experimentally achievable parameters. 
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II. FUNDAMENTAL LIMIT FOR PHOTON 
BLOCKADE WITH A TWO-LEVEL SYSTEM 

The first demonstration of the photon blockade was 
performed with Cs atoms strongly coupled to a Fabry- 
Perot cavity [3], and later with quantum dots coupled 
to a photonic crystal resonator [6j. In this strong cou- 
pling regime, where the coupling rate between the atom 
and the cavity is larger than all the other rates in the 
system, the photon blockade relies on the anharmonic- 
ity of the energy levels of the coupled atom-cavity sys- 
tem [U [TO], meaning that multi-photon excitations of 
the coupled system occur at different frequencies than 
a single-photon excitation. A similar anharmonicity of 
the spectrum of semiconductor systems has been used to 
demonstrate single-photon sources [IHl HI] ■ 

The second photon blockade mechanism occurs at the 
bad cavity limit [22', in which the cavity-enhanced cou- 
pling of the atom to one electromagnetic mode is faster 
than its spontaneous emission to all other modes, yet 
slower than the cavity decay rate. Therefore, although 
strong coupling is not achieved, the atom interacts mostly 
with one electromagnetic mode and can thus be perceived 
as a one-dimensional atom [53]. The blockade effect is 
then typically described as a dynamical process, equiv- 
alent to antibunching in free space resonance fluores- 
cence IM' . This stems from the fact that in the bad cavity 
regime, the scattered photon escapes the cavity immedi- 
ately, before the atom, which collapsed to the ground 
state, can rebuild its polarization and scatter another 
photon. Nonetheless, one can consider this mechanism in 
the spectral domain as well, like in the strong coupling 
regime. The difference is that in this case the anhar- 
monicity is of the atomic energy levels, where, again, 
multi-photon excitations occur at different frequencies 
than a single-photon excitation. In either case, the fi- 
nal result can be described as a "photon turnstile" in 
continuous operation, since photons are transmitted (or 
reflected) one by one, as was demonstrated with Cs atoms 
coupled to a whispering gallery mode of a fiber-coupled 
microtoroid cavity (S] [11] . 

A. Photon blockade in the bad cavity limit 

We begin by analyzing the blockade mechanism in the 
bad cavity limit. For this we consider a system similar to 
that of Ref. [TT], namely a single atom interacting with 
a fiber-coupled microtoroid, as depicted in Fig. [T] 

The two counterpropagating modes of the microtoroid, 
d and b, are evanescently coupled to the atom with cou- 
pling strength g, and decay to the tapered fiber at a rate 
2Ke^ 3> g. Since we are interested in inherent limits to 
the routing effect, we shall assume at first that parasitic 
effects such as scattering between the cavity modes, 
intrinsic losses of the cavity and the natural decay of the 
atom are significantly slower than all the other processes 
in the system. The dominant parameter in the system is 



thus the cavity decay rate into the fiber, which brings us 
into the bad cavity regime. It is important, however, to 
make a distinction between this situation, in which the 
cavity is intentionally made to decay into the fiber, and a 
truly "bad" cavity that decays quickly due high intrinsic 
losses. Thus, we shall henceforth refer to this regime as 
the "fast cavity" regime, instead of the bad cavity regime. 

We define the routing efficiency C*'" as the probabil- 
ity of detecting one photon in the transmitted mode ciout 
and one in the refiected mode bout, given that the input 
mode din initially contains a two-photon pulse. Neither 
the time interval between the two clicks, nor their ab- 
solute time of occurrence are of importance, and hence, 
for calculating the probabilities for the four possible de- 
tection events, we integrate over all possible detection 
times; 
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where r stands for refiected and t for transmitted, and 
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are the second order correlation functions of the two out- 
put modes. Note that this definition is already normal- 
ized so that P*'' + P''* + P''' + P" = 1. The appropriate 
expression for the routing efficiency is thus given by 

Accordingly, P** and P'"'" denote the probabilities of the 
two failure mechanisms of the routing process, namely 




FIG. 1: (Color online) Schematic depiction of a microtoroid 
setup. Note that bin (not shown here) is vacuum. 
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FIG. 2: A schematic of the photon router model. The input 
pulse is modeled by introducing a feeder cavity which leaks 
into the system with rate 2ks. 



The non-Hermitian effective Hamiltonian correspond- 
ing to the cascaded system of Fig. [2] is given by {h = 

1) m 

Hi = —iKgiilag — 2^7^(7^(7 — 2i^Ks7c(7^a<,. (5) 

Note that the interaction between the feeder cavity and 
the system is unidirectional, and that it produces an 
exponentially decaying driving pulse. Also, since the 
atom can decay into two counterpropagating modes, the 
atomic population decay rate is instead of 2^^- 



C. Results 



the probabilities for both photons to be transmitted, or 
for both photons to be reflected, respectively. 

For an ideal router, events where both photons are 
transmitted or both reflected should not occur: P** = 
prr _ g^j^j ^]-^^g (jtr _ ^ comparisou, for a sim- 
ple 50 : 50 beam splitter, any of the four possibilities is 
equally likely, hence C*'' = 0.5. Any useful photon router 
should thus satisfy C*'" > 0.5, which is also the quantum 
limit, above which the Cauchy-Schwartz inequality for 
classical fields is violated. 



B. One- and two-photon source: the feeder cavity 

In order to analyze the temporal behavior of a pulse 
interacting with this system and the resulting routing ef- 
ficiency, we model the fiber by a feeder cavity containing 
the desired number of excitations (Fig. [2]). 

The decay rate of the feeder cavity 2ks gives a measure 
of the inverse of the resulting driving pulse width. Using 
the input-output formalism |25H27| and eliminating the 
microtoroid cavity modes adiabatically, which is justified 
in the fast cavity regime for which Kex 3> g, we can write 
for the output modes 



^out 

bout 




(4a) 
(4b) 



where ds is the annihilation operator of the feeder cav- 
ity field, a is the lowering operator for the atom, and 
27c = 2g'^ /Kex is the cavity-enhanced atomic decay rate 
per output channel. Since only normally-ordered mo- 
ments of the output operators are considered, the vacuum 
noise input operators fls.m and bi„ have been discarded. 



Eq. (4a) reflects the fact that emission of a photon into 



&out by the feeder cavity and by the cavity-enhanced 
atom are indistinguishable. Accordingly, it is the interfer- 
ence between both paths that is detected. However, since 
we assumed the scattering between the cavity modes is 
negligible, only the atom is able to reflect a photon, and 
hence no interference will be observed in the reflected 
mode described in Eq. ( 4b I (assuming is vacuum) . 



1. Single photon driving pulse 

As a simple example, consider a driving pulse contain- 
ing a single photon. The initial state is then 



l^(0)> = |lff> , 



(6) 



where the first index in the ket gives the number of pho- 
tons in the feeder cavity, and the second index describes 
whether the atom is excited (e) or in its ground state (g). 
The state evolves according to the Schrodinger equation 
iijj = Hiip to 



\m)^a{t)\lg)+b{t)\Oe), 



with 



a{t) 
b{t) 



(7) 



(8) 



2y/Kslc 

27c - 



) 



where henceforth the exponential functions are defined 
to be zero for t < 0. We use Eq. Q to evaluate the 
photon flux at the output modes with this state, and 
integrate over time to define single-photon transmission 
and reflection probabilities 
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- (9a) 
, (9b) 



respectively. 



As evident from Eq. (|9]), for long enough pulses (k^ ^ 
7c), T — T' 0, i? — T' 1 and all light is reflected. This is due 
to the destructive interference between the atomic dipole 
radiation and the driving field in the forward direction. 
For short pulses (ks ^ 7c), however, the atom has no 
time to build up the necessary dipole, and the photon is 
transmitted. 
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2. Two photon driving pulse 

For deducing the routing properties of this system, we 
need to study the deviation from the hnear behavior of 
Eq. (|9| and analyze the case of more than one photon 
in the pulse. The source of the nonlinearity is evident 
already in Eq. Q, in the presence of the atomic low- 
ering operator in the expressions for the output mode 
operators. In particular, detection of a reflected photon 
projects the atom to its ground state, preventing the im- 
mediate scattering of a second photon and also disrupt- 
ing the destructive interference in the forward direction, 
possibly allowing the second photon to slip through and 
be transmitted |28j . In order to derive an analytic de- 
scription for this effect we follow the evolution of a state 
in which the feeder cavity initially contains two photons. 
This initial state evolves to 

= 125) +/? Wile), (10) 



where 

a[t) = e-2«=t (11) 
^^^^ ^ ^ y27c^s [^^..t _ ^-2^.t] ^-.^t_ 

After the detection of one photon, the wave function col- 
lapses to a state containing only a single excitation, which 
can be either in the cavity or in the atom. In the first 
case the subsequent evolution is again given by Eqs. Q 
and ([S]), whereas in the second case, the state evolves to 
|^(t)) = c{t) |0e), with 

c{t) = e-27ct, (12) 

Substituting these expressions into Eq. ([2|, the second 
order correlation functions become 
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(13) 



Substituting these correlation functions into Eq. ([T]) and 
performing the integration leads to the following expres- 
sions for the various probabilities for transmission or re- 
flection of the photons 
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These probabilities and the corresponding routing efh- 
ciency C*'' = P*^ + P^* are presented in Fig. [s] as a func- 
tion of the pulse width. 

As evident, there is no pulse width for which both fail- 
ure mechanisms P**, P^^ are zero, and so the routing ef- 
ficiency C*'' is a compromise, reaching a maximal value 
of only ~ 64%. These results are consistent with pre- 
vious works by Koshino et al. [3D] and Shapiro [3T] , 
although our model and derivation are different. The 
physical origin of the limited routing efficiency relies on 
the interplay between two counteracting effects. On one 
hand, if the pulse is significantly longer than 7~^, the 



atom has time to reestablish its dipole after the first 
scattering, and therefore will be able to scatter the sec- 
ond photon backward as well when it arrives. On the 
other hand, due to time-energy uncertainty relations, a 
pulse that is significantly shorter than 7"^ must have 
a bandwidth that exceeds 27c, which is also the interac- 
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FIG. 3: (Color online) Two-photon detection probabilities as 
a function of the pulse width kJ^. 
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FIG. 4: (Color online) Marginal probability density functions 
of the various reflection and transmission events, as a func- 
tion of the interval between the photon detection times r, for; 
(a) long pulses {ks — 0.057c), (b) short pulses (ks = 57c) and 
(c) intermediate pulses (ks = 1.57c). The different curves 
correspond to events of double reflection (RR), double trans- 
mission (TT), a transmission followed by a reflection (TR), 
and a reflection followed by a transmission (RT). Note that 
RR is practically zero in (b). 



tion bandwidth of the atom. Thus, such a pulse will have 
some spectral components that cannot interact with the 
atom, and consequently there will be a probability for 
both photons to be transmitted. Since these two effects 
occur on the same timescale, the efficiency is inherently 
limited. 

In order to gain more insight into the dynamics of this 
process, we present in Fig. |4] the correlation functions of 
Eq. ( 13 1 with t integrated over the entire pulse for inter- 
mediate, long and short pulses. This integration results 



in a marginal probability density function, displaying 
the relative probability for each of the four detection 
events as a function of the time interval between the two 
detections r. A second integration over r would result 
in the probabilities of Eq. ( 14 ) . 



For all pulse lengths the reflection-reflection compo- 
nent RR is antibunched, corresponding to the projection 
of the atom to the ground state following the reflection 
of the first photon. This is the blockade effect, which 
implies that photons are reflected on a one-by-one ba- 
sis. The antibunching of RR naturally leads to a com- 
plementary peak in the reflection-transmission graph RT, 
since at that timescale a second photon must be transmit- 
ted. The transmission-transmission component TT, how- 
ever, results from different mechanisms at short and long 
pulses; at short pulses it reflects the probability that the 
(spectrally broad) photons did not interact with the atom 
at all, and so the peak in TT is as narrow as the pulse. 
For long pulses, TT results from the probability that the 
atom let one photon 'slip through', since it was already 
excited at that moment by the other photon. In that case 
the atom can emit the second photon both forward and 
backward at equal probabilities; emission forward is the 
one that is presented in the TT graph, which is thus as 
wide as the atomic time scale 7~^. Backward emission 
of the second photon leads to a transmission-reflection 
event TR, which is perhaps the most interesting process 
of the four. In particular, as we see in Fig. |4j;, TR is the 
component that contributes the most to the routing effi- 
ciency C*'' at intermediate pulses, and not, as one would 
have expected naively, the antibunching of the RR pro- 
cess, whose contribution is negligible at intermediate and 
short pulse lengths. 

The underlying mechanism of the TR event with inter- 
mediate pulses is the simple fact that the photon emis- 
sion rate from a feeder cavity occupied by two photons 
is twice as large as the emission rate from a cavity con- 
taining only one photon. Consequently, a first photon 
is very likely to arrive shortly after the beginning of the 
pulse, before the atom builds its dipole field, and thus 
it is likely to be transmitted. The time interval between 
the first photon and the arrival of the remaining second 
photon is expected to be longer, giving the atom enough 
time to reflect it. As in the case of TT, transmission 
of the first photon does not necessarily leave the second 
photon in the feeder cavity, but can also result in the 
collapse of the atom to its excited state, from which the 
photon has a 50% probability of being refiected. These 
two mechanisms interfere constructively for the case of 
TR, but destructively for TT, thus considerably enhanc- 
ing the routing efficiency. 



D. Simulations 

In our calculations so far, the temporal profile of the 
pulse is the only free parameter in the problem, and thus 
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the limited efficiency of the photon blockade mechanism 
is inherent and unavoidable. In order to verify the an- 
alytic calculations, and also to check the effect of using 
various pulse shapes, a wave function approach using a 
nonunitary Hamiltonian was used to provide a fully quan- 
tized and complete simulation of the setup, including the 
microtoroid cavity (which was adiabatically eliminated 
in our analytical calculations) . The initial two-photon 
state is specified by 



al{xi)al{x2) 
V2 



\0g) dxidx2, (15) 



where a\.{x) and a^{x) {Sl[{x) and a;(a;)) create and anni- 
hilate a right (left) propagating photon at a location x in 
the fiber. |0) denotes the vacuum state of the microtoroid 
cavity, and \g) the ground state of the atom. f{xi,X2) 
is a normalized weight function that describes the prob- 
ability amplitude of the two photons to be located at xi 
and X2 ■ The effective Hamiltonian is given by [32l [33] 



dxdl{x) ^wo — ar{x) (16) 

-I- j dxa]{x) {^o + iVg^^ ai{x) 

+ J dxxix) {Vdl{x)a + V*a^ar{x)) 
+ J dxxix) (Va]{x)b + V*b''di{x)'^ 

with d and d^ (b and W) the annihilation and creation 
operators associated with the counterclockwise (clock- 
wise) microtoroid cavity modes of frequency ujc- is 
the central frequency of the pulse around which the fiber 
dispersion relation is linearized. a\ and ag are the bare 
atom raising and lowering operators of the transition with 
frequency tUa ■ Vg is the group velocity of the pulse as it 
propagates through the fiber. 27 is the population decay 
rate of the bare atom and 2/*^ is the intrinsic microtoroid 
cavity decay rate. The coupling strength between the 
fiber and the cavity is denoted hy V = y/2KexVg. The 
coupling between the cavity and the fiber is assumed to 
have a normal distribution xix) cx e~2(^^/^'^'> , where 
Lt is the effective interaction length. 

Figure [5] presents the routing efficiency for various in- 
put pulse shapes, demonstrating little dependence on the 
pulse profile. Gaussian pulses have been found to yield 
the best results, with a routing efficiency of C*'" — 66.8%, 
yet the inherent conflict and the resulting limited effi- 
ciency remain. 
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FIG. 5: (Color online) Simulated results for the two-photon 
detection probabilities as a function of the pulse width for 
square pulses (dashed curves), Gaussian pulses (dash-dotted 
curves) and exponential pulses (solid curves). For square and 
Gaussian pulses the pulse width k,~^ is defined by their half 
width half maximum. The parameters are g/2n = 70 MHz, 
Kexf^n =500 MHz, without losses (k; =7 = 0). We note that 
the simulated results for exponential pulses are practically 
identical to the analytically calculated results presented in 
Fig-i 



E. Photon blockade in the strong coupling regime 

The analysis so far relied on the fact that only 
one-photon transitions could take place, since atomic 
two-photon transitions (such as from 55*1/2 to 5D^/2 
in Rb) occur at a different frequency. Thus, the atom 
behaves as a two-level system capable of dealing with 
only one excitation, with a linewidth dictated by the 
cavity-enhanced coupling to the output modes 7c. 
Exactly the same results are expected in the case of 
photon blockade in the strong-coupling regime of cavity 
QED. In that regime, the nonlinear system is not the 
atom, but rather the coupled atom-cavity system, in 
which one-photon resonance and two-photon resonance 
are indeed at different frequencies [31]. The linewidth 
of each of the vacuum Rabi sidebands in that case is 
[k + 7)/2, with K being the bare cavity linewidth, and 
7 the atomic linewidth in free-space. Thus, all our 
results so far are relevant for the photon blockade in the 
strong coupling regime as well, in which the interaction 
involves one of the Rabi sidebands, at detuning g from 
the bare-cavity resonance. The only difference therefore 
is that the bandwidth (k + 7) /2 replaces 7c in all the 
preceding expressions. We verified this by performing 
fully quantized simulations using the above method in 
the strong coupling regime, with the photons tuned to 
one of the Rabi sidebands. Indeed, these simulations 
yield results indistinguishable from those obtained in 
the fast cavity regime. 
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FIG. 6: The feeder cavity is driven by a three-level atom 
emitting two photons at once. 



III. ROUTING EFFICIENCY USING 
TIME-ENERGY ENTANGLED PHOTONS 

As established in the previous section, time-energy un- 
certainty relations are the key mechanism that limits the 
routing efhciency. Therefore, it is interesting to study the 
possible influence of using time-energy entangled photon 
pairs as the input two-photon pulse. 



A. Three-level atom source of photon pairs 

In order to obtain analytic expressions, the input pulse 
is now modeled by introducing a two-photon emitting 
atom with an infinitesimal intermediate state lifetime 



(Fig. [6]) driving the feeder cavity. Two photons are thus 
simultaneously emitted into the feeder cavity at a rate 
27s, from which they decay independently at a rate 2k, g. 
The feeder cavity then 'smears out' the ideal entangle- 
ment generated by the three-level atom. Therefore, the 
ratio Kg lis sets the amount by which the two-photon 
pulse that drives the router is temporally squeezed. 
The effective Hamiltonian to be used in this case is 

H = —iKga\ag — 2i^s(jg^ dg — 2«7c(7^<t (17) 
-2iy/2K'gjg(fg (alY - 2iy^Kg-fcd-''ag, 

where ds\ dg are the raising and lowering operators of 
the three-level atom for the transition from the upper 
state e to the ground state g. n'g <^ Kg is the linewidth 
of the left mirror of the feeder cavity, and is taken to be 
very small but nonzero, to allow the feeder cavity to be 
driven by the three-level atom. Note that although the 
quantity of interest is the entangled pulse transmitted to 
the router, most of the pulse does not enter the feeder 
cavity, and hence only events whereby both photons enter 
the feeder cavity are taken into account. The two-photon 
state can be expanded as 

m)) = m |eOg) + a{t) \g2g) + /3(i) \gle) , (18) 

where the first index describes the state of the three-level 
atom. Solving the Schrodinger equation we obtain 
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(27, - Kg)e-^''^' + (27, - 27, - Kg)e-^^^' + 2{Kg - 7,)e-(2T^+"^)' 
(27c - Kg){2jg - 27c - Kg){Kg - 7s) 
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By substituting these expressions into Eq. ( 13 1, the rout- 
ing efficiency can be calculated using Eqs. (l]|3l, where 
a(t), b{t) and c(t) are those of Eqs. ( 8|12 1. This leads to 
the following expressions for the transmission and reflec- 
tion probabilities: 
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(20) 



Note that these equations coincide with Eqs. (14) in the 
limit of an unentangled input pulse (7^ — 00). We see 
that the routing efficiency (Fig. [7| increases significantly 
as the entanglement is increased, reaching a 10% im- 
provement at Kg/^g « 5. For even larger values of Kg/^g, 
the increase becomes less pronounced, and the efficiency 
tends to an asymptotic value of ^ 77%. Thus, we see that 
even an entangled state at the input cannot completely 
eliminate the conflict that limits the routing efficiency of 
a two-level system. 



B. Cavity QED source of photon pairs 

Although an entangled two-photon input pulse does 
not enable ideal routing, the increased efficiency obtained 
in this scenario seems potentially useful, especially since 
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0.5 1 1.5 2 

(units of y~^) 

FIG. 7: (Color online) The routing efficiency C*"^ as a function 
of the pulse length k.J^ and the ratio Ks/Js- 



time-energy entanglement can be generated within a two- 
photon pulse by reflecting it from a single-sided cavity 
containing a two-level atom [T5l[TS], as depicted in Fig.|8] 
Thus, one could imagine performing a two-stage inter- 
action with a cavity QED system, the first for gener- 
ating the entanglement, and the second for performing 
the routing. The underlying physical effect behind the 
generation of entanglement is the same interference that 
occurs between the driving field and the atomic dipole, 
which suppresses the forward transmission in the regu- 
lar cavity configuration, and leads to the backscattering 
of the photon. However, in the case of a single-sided 
configuration, there is only one cavity mode, and the de- 
structive interference between the atomic dipole and the 
driving field is only partial, such that the emitted photon 
flux equals the incoming photon flux. Once a photon is 
scattered, however, the atomic state flips its phase, and 



[Feeder cavity] 



QWP 



PBS 



7c2 



Jell 




^out 



FIG. 8: Cascading of a single-sided cavity QED system and 
a photon router. Due to the quarter wave plate (QWP), the 
pulse is reflected by the polarizing beam splitter (PBS) after 
interaction with the single-sided cavity. The trajectory of the 
two-photon pulse is shown by the dashed line. 



subsequently the destructive interference becomes a con- 
structive one. In more detail, the effective Hamiltonian 
modeling the passage of a pulse through a single-sided 
cavity QED system in the fast cavity regime is given by 



(21) 



where we assume that the atom is coupled to only one 
standing wave in the cavity, explaining why the decay 
rate of the atom is twice as small as in Eq. ([5|. The 
two-photon state is now expanded as 



|^(t))=a(t)|25)-t-/3(i)|le). 



with 



a{t) -- 



e-7citl g-Kst_ 



7ci 



(22) 



(23) 



Note that due to the single-sided conflguration of the cav- 
ity, these expressions differ from those of Eq. (111. It is 



evident that the driving and the atomic radiation are out 
of phase. To illustrate the dynamics more clearly let us 
assume a long pulse (ks <^ 7ci) and neglect transient ef- 
fects that occur at t ^ 7^j^^ . Under these approximations 
the state of the system is 



\m) 



\2g) -2V2J—\le) 
7ci 




Applying the output operator of Eq. ( 4a ) yields 



aout Hit)) « 



(24) 



(25) 



corresponding to a photon detection probability rate of 



4k., 



-4K.st 



As evident, the detection of 



the first photon has led to a sign flip between the driving 
field and the atom, resulting now in a constructive 
interference, rather than the previous destructive one. 
Normalizing the state and applying aout again shows 
that this leads to a sudden increase of the photon detec- 
tion rate by a factor of 4.5 to ^ a^^^hout^ ~ 18 Kse~^'^'* . 
Since after one photon detection there is only one 
excitation left in the system, naively one could have 
expected the probability for a second photon detection 
to drop and not increase; this increase is therefore in 
fact by a factor of 9 compared to the expected value at 
steady state, and thus the output state exhibits strong 
bunching and consequently time-energy entanglement 
between the two photons [l6] . 



Solving the evolution of the two-photon state ( 22 ) us- 



ing Eq. (21) without approximations, we get that the 



probability amplitude of photon detections at times t and 
t + Tis (Fig. l9| 



fit.r) 



1 - 4e-('''^i"«='^ 



(26) 



-2 1 - 3e"(^'^i"'"='^ e-^citg-«,tg-K,T 
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FIG. 9: (Color online) Probability of detecting two photons 
at times t' and t" for an exponential pulse {"fci — 5ks) after 
passage through a single-sided cavity QED system. 
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y^2 (units of K^) 



FIG. 10: (Color online) Routing efhciency as a function of 
the cavity-enhanced decay rates during the first and second 
passages through the system. 



The terms in the second brackets correspond to the 
transient behavior of the cavity-enhanced atom, and 
have influence only for short pulses, as evident from 
their fast decay by e~'^'''^*. The two terms in the first 
brackets thus demonstrate the main dynamics of the 
system. The first term presents the possibility that no 
interaction occurred with the atom, and the second, 
nonlinear term is the one that leads to the bunching of 
the two-photon pulse. 

The entanglement appearing in Fig. [9] indeed suggests 
that by utilizing this pulse as the input to a routing 
system, higher efhciency might be obtained. In order 
to investigate this possibility, a triply cascaded system 
needs to be considered (Fig. [s]). The feeder cavity 
drives a single-sided cavity QED system with effective 
linewidth 7ci, and the resulting bunched light drives the 
photon router with linewidth jc2- 



The Hamilonian of this system consists of the Hamil- 
tonians of its subsystems, given by Eq. ^ and Eq. 
and a term corresponding to the driving of the photon 
router by the single-sided cavity 



H3 ~ Hi + H2 



+ VtcK^i). (27) 



The analytic solution of the detection probabilities 
is calculated in the appendix. The resulting routing 
efficiency is shown in Fig. |10[ 

Note that for jci = we recover the result of Fig. [3] By 
increasing jci the pulse becomes bunched, and the rout- 
ing efficiency reaches ^ 68%. However, while the bunch- 
ing ratio increases, the bunching efficiency decreases. 
This is due to the decreasing area of the bunched part 
in Fig. |9] as it becomes narrower. Hence, further in- 
crease of 7ci ceases to improve the routing efficiency. We 
conclude that a two-stage interaction with the cavity- 



enhanced atom, can only improve the routing efficiency 
by at most ^ 4%. 



IV. IDEAL ROUTING EFFICIENCY USING A 
THREE-LEVEL ATOM 

In the previous sections we have studied the limitations 
of a two-level system as a photon router. The situa- 
tion changes dramatically when another cavity-enhanced 
atomic transition is introduced, for example by utilizing 
a three-level atom in the A-configuration inside the cav- 
ity. We shall consider the behavior of this system in the 



single-sided cavity setup, as depicted in Fig. 11 

The two ground states are denoted by \gi) and I52), 
and the excited state by |e). We denote the fields res- 
onant with the gi ^ e transition and with the 52 e 
transition as H and V fields, respectively, assuming they 
differ in polarization. This enables separation of both 
fields by a polarizing beam-splitter (PBS). In this con- 
figuration, therefore, routing to the two different output 
ports is performed by manipulating the polarization of 
the photons 

Taking the feeder cavity to be resonant with the H- 
transition, the appropriate Hamiltonian is 



PBS 




ay 




FIG. 11: Schematic depiction of a three-level atom in the A- 
configuration inside a single-sided cavity. The trajectory of 
the pulse is shown by the dashed line. 
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(28) 

where aee = |e)(e| and aie = |5i)(e|- '2^cH, 2jcV are 
the cavity-enhanced decay rates from the excited state 
into ground states \gi) and \g2), respectively. 

For a long input pulse (ks <^ IcH^lcv) containing a 
single i?-photon, the probability of a full transfer to the 
right ground state \g2) starting with an atom in the left 
ground state \gi) is given by 



(29) 



1 + {1cV-1ch) I'^lcVlcH 



Thus, by choosing transitions for which ^cH — IcV = 
one can approach Py — !■ 1, and obtain determinis- 
tic transfer of the atom from one ground state to the 
other. In that process, the iJ-photon is absorbed and 
released as a V^-photon. This configuration was studied 
by Koshino et al., who demonstrated that it may be used 



to implement deterministic quantum state transfer and 
a VSWAP gate between a photon and an atom [17l |35] . 

In order to use this apparatus for photon routing, we 
generalize this result by starting with a Fock state |n) 
of type H in the feeder cavity. Assuming the atom is 
initialized in the left ground state \gi), the wave function 
evolves to 



m)) =a{t)\n,gi)+l3{t)\n~l,e) 



(30) 



By solving the Schrodinger equation with the Hamilto- 
nian of Eq. (28), we get 



a{t) = 6-""=* 



(31) 



m = 



IcH + IcV - 



Using the output operators from Eq. Q , we get that the 
probability of first detecting a ^-photon is 



Pv{n) 



1 

{icV - Icnf + {icv + Ich) (3rt - 2) + {2n - 1) (n - 1) nl j^^icHlcV 

I 



(32) 



Once again, by using a long enough pulse and sym- 
metric transitions, the first detected photon will with cer- 
tainty be a l^-photon. When that happens, the atom col- 
lapses to the ground state \g^l)^ and becomes transparent 
to all the remaining i7-photons, which are consequently 
reflected from the cavity unchanged. The mapping re- 
sulting from the passage through the system is thus 

(aL)"|0)^ (4)""'«y|0>, (33) 

i.e., it acts as an ideal photon router, which is transpar- 
ent to all but one photon. 

Specifically, in the case of n = 2, we obtain for the 



routing efficiency (Fig. 12 1 



HV, 



(34) 



where Pvh = Pvin — 2) is the probability of detect- 
ing an iJ-photon after a ^-photon has been detected, 
and Phv is the probability of first detecting an H- 
photon and then a V^-photon. Taking into account that 
Pvh — Pv — 2) can be made close to unity by itself, we 
see that this routing process includes an inherent robust- 
ness, or error-correction mechanism, since in the unlikely 
event in which the first photon stays H (for example, due 
to a slight mismatch between JcH and jcV or due to non- 
negligible Ks), the second photon is likely to be turned 
into V, leading to the desired result nonetheless. Naively, 
one could expect Phv to be {1 — PvH)Pv{n =1), namely 



the probability that the first photon was not turned to V 
times the probability of turning the last remaining pho- 
ton into however Pv{n) is correct only for an atom 
that begins at the ground state \gi), whereas after de- 
tection of an H photon the atom has some probability 
amplitude to be in the excited state |e). The full solu- 
tion for Phv, following the same procedure as for Pvh, 




FIG. 12: (Color online) Routing efficiency C*"^ as a function 
of the cavity-enhanced decay rates "/cH and "/cv- 
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Phv = 7 , , ^7 , ^ Pvin = 2). 

(35) 

The overall routing efficiency is presented in Fig. [12] 
The contribution of Phv adds to already-high values 
of PvH to result in efficiencies that approach unity for 
reasonable jcV^lcH, and thus we see that this process 
acts as a robust and nearly ideal photon router, free 
from the limitations discussed in the previous sections. 

Finally, we wish to estimate the robustness of this rout- 
ing scheme to losses in the system. We approximate the 
effect of the coupling to the environment by calculat- 
ing the probability of a quantum jump in the system, 
neglecting the change in Eq. (31). The effect of sponta- 



neous emission into the vacuum at rate 7 for long pulses 
and jcH = IcV = 7c is then given by 



psp 



(36) 



which is independent of k^. The efficiency of the routing 
is thus not increased by taking longer and longer pulses, 
as is the case for stimulated Raman adiabatic passage, 
but can be suppressed by increasing the cavity-enhanced 
atomic decay rate. The effect of intrinsic losses of the 
microtoroid cavity can be evaluated by omitting the fast- 
cavity approximation, and including the microtoroid cav- 
ity modes in the system dynamics. The probability of 
losing a photon is then calculated to be 



(37) 



This effect can be suppressed by increasing the cavity 
decay rate into the fiber. Hence, by taking 7c ^ 7 and 
Kex 3> Ki, the routing mechanism is made robust against 
both loss mechanisms. 



system is insufficient for this task, since the bandwidth 
of the system and the timescale of its nonlinearity, 
namely the temporal duration at which the system will 
respond differently to two photons than to one photon, 
are inherently linked by the uncertainty principle. 
This inherent conflict limits the efficiency of photon 
routing to ~ 64%. We have shown that even the use of 
time-energy entangled photons at the input pulse does 
not circumvent this conflict completely, although it does 
enable some increase in the efficiency of the process. 
Finally, we have shown how the use of a three-level 
system does enable ideal photon routing, free from the 
limitations of a two-level system. Specifically, we pre- 
sented a scheme that uses an atom in the A-configuration 
with both transitions enhanced by a single-sided cavity, 
to create a photon routing mechanism in which one and 
only one photon from the incoming pulse is directed to 
one port, while the remaining photons are directed to 
the other port. This scheme is robust against variations 
in the pulse width and parasitic losses with realistic pa- 
rameters, and provides a promising method for efficient 
routing of single photons. 
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Appendix A 

The two-photon state of the cascaded system depicted 
in Fig. [8] can be decomposed into 



V. CONCLUSION 

In this work we analyzed photon routing in cavity- 
QED, focusing on the task of separating two incoming 
indistinguishable photons to different ports. We have 
shown that using just one cavity-enhanced two-level 



\m) = |2g.9) + m I leg) + Sit) \lge) + r^{t) |Oee) , 

(Al) 

where the second and third indices describe the atomic 
states of the single-sided subsystem and the photon 
router subsystem, respectively. Using the Schrodinger 
equation, the coefficients of Eq. ( Al I are found to be 
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a{t) 

m 

Sit) 

vit) 



-2k. t 



7cl - '- 



(A2) 



„-(27c2+Ks)t 

7cl - Ks 27c2 - Ks 7cl - Ks 2jc2 ~ 7cl 27c2 - 7cl 27c2 ~ 

27c2 - 7ci - 2ks 



4^2 



«^sV7cl7c2 
(7cl - Ks)(27c2 - Ks) { 27c2 + 7ci - 2k 



g-2K,i _ g-(27e2+7cl)t 



37cl - 27e2 

27c2 - 7c 1 



g-(7cl-l-Ka)t _ g-(27c2+7cl)t 



27c2 + 7ci 
27c2 - 7ci 



g-(27c2+Ka)i _ g-(27c2+7cl)t 



After a photon detection, the system can collapse into 
three different single-excitation states. If the photon col- 
lapses to the feeder cavity, the state evolves according 
to 



\Mt)) = ai |l53> + bi \Oeg) + ci \Oge) 



where 



(A3) 



where 



b2{t) 
C2W 



2 V7cl7c2 . 

27c2 - 7ci 



-27c2n 



(A6) 



aiit) 

hit) 

ciit) 



-2 (e-«=* 

7cl - 



^ V7c2Ks j 7cl +K^ ^-K,t 



27c2 - Ks [ 7cl - K 
27c2 + 7cl ^-27c2t 

27c2 - 7ci 



(A4) 



2 c~^"'* 
7ci - 



If the system collapses to a state in which the first 
atom is excited, we have 



\Mt))=b2\0e9)+C2\09e), 



And finally, if the state collapses to a state in which the 
second atom is excited, we have |'03(^)) = c^it) \Oge), 
with 



C3(t) = 



(A7) 



(A5) Using Eq. Q, the correlation functions are given by 



r*'-(t,T) = 4\Vl^vit)Wlc.lb2{T)+^C2{T)]+^S{t)[^sai{T)+y/^bl{T) + y/^Cl{T)]f (A8) 

r'-'it, t) = 4 I [V2^«(t) -I- V^/3(i) + Vl^6{t)] ci(r) 

+\/7^ [V7^77(i) + \/Kr,5(i)] C2(t) + \/>2 [\/>2'7(*) + %/K75(i)] C3(t)|^ 
T^^{t,r) = 4\j,Av{t)c2{T) + S{t)c^{r))f 

r"(i, r) - 4 I [x/2^a(t) + + V7^(5(t)] [y^ai(r) + V7^&i(r) + V7^ci(r)] 

+ [Vl^v(t) + v^sPit)] [V7^62(t) + V7^C2(t)] -I- 7^ [V7^?7(t) + V'i^sS{t)] C3(r)|^ . 



Integrating these expressions yields the results plotted in 
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